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ABSTRACT 

Observations of the globular cluster 47 Tucanae (NGC 104), which contains at least 23 millisecond pulsars, were performed with the H.E.S.S. 
telescope system. The observations lead to an upper limit of F(E > 800 GeV) < 6.7 x 10"'^^ cm"^ s"' on the integral y-ray photon flux from 
47 Tucanae. Considering millisecond pulsars as the unique potential source of y-rays in the globular cluster, constraints based on emission models 
are derived: on the magnetic field in the average pulsar nebula and on the conversion efficiency of spin-down power to y-ray photons or to 
relativistic leptons. 

Key words. (Galaxy:) globular clusters: individual: 47 Tucanae - (Stars:) pulsars: general - Gamma rays: observations - msPSR 



1. Introduction 

Millisecond pulsars (msPSRs) are usually categorized among the radio 
pulsar population by limits on their spin period (P < 50 ms) and, when 
available, intrinsic spin-down rate (P,,,, < 10""* s s"'). They are old neu- 
tron stars, possibl y re-acceler a ted by interactions with a companion, as 
first proposed in lAlpar et al.l ( Il982h . Very high energy (VHE) emis- 
sion from this type of object has been predicted via various radiation 
mechanisms. For individual objects. Inverse Compton (IC) or Curvature 
Radiation (CR) emission due t o the acceleration of leptons above the 
polar cap jHarding et al.ll2005l : iBulik et ani2000h have been proposed. 
For binary systems, an additional possibility would be the interaction 
between pulsar wind driven outflows and the stellar wind of the com- 
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panion (see for instance rDubusll2006h . The spin-down power, typically 
lower than 10^^ erg s"', entails expected individual y-ray fluxes well be- 
low the detection threshold of current instruments. 

However, groups of msPSRs have been identified in Galactic glob- 
ular clusters (see e.g. lManchester et aljl99ll) . allowing for larger fluxes 
from an ensemble of unresolved sou rces. Out of the 185 pu lsars with 
P < 50ms known in the year 2008' ^Manchester et al.ll2005l) . 140 be- 
long to globular clusters^. Globular clusters (GCs) are old high-density 
galactic str uctures, with ages c lose to the age of the Galaxy itself (see 
for instance fGratton et al.l2003l) . Their age indeed suggests evolved em- 
bedded stellar populations including compact (binary) objects, which 
are considered potential progenitors to msPSRs, as discussed e.g. in 



' http://www.atnf.csiro.au/research/pulsar/psrcat/ [vl.34] 

- http://www2.naic.edu/~pfreire/GCpsr.html [on 2008 August 7] 
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iBenaguistal ( 120061) . The GCs Terzan 5, 47 Tucanae and M28, in this 
order, host the largest identified msPSR populations teanso n]|2008h . 

47 Tucanae (NGC 104) is one of the largest Galactic GCs known 
to date, with an estimated mass of 10' Mq and an age of 11.2 + 
1.1 Gyr jGratton et alj2003|). Optical observation s by the Hubble Space 
Telescope, described in iMcLaughlin et alj ( l2006h . allowed precise mea- 
surements of its location, centered at aaooo = 0'' 24"' 05 '.67 and <52ooo = 
-72° 04' 52".62 and placed it at a distance of 4.0 ± 0.35 kpc. The sur- 
face brightness distribution allows the estimation of a core radius of 
ro = 20". 84 ± 5". 05, a half- mass radius r;, x 2.6' and a tidal radius 
r, » 0.6°, using the model of iKind ( 119661) . In 47 Tucanae, 23 pulsars 
so far were revea led, with radio ob servations predominantly using the 
Parkes telescope dFreire et al.l2003l) . with periods in the range 2-8 ms, 
averaging at 4 ms, all located within 1.2' of the centre of the GC. 
Based on the unresolved 20 cm radio flux from the core of 47 Tucanae, 
iMcConnell et alj ( 120041) estimated that up to 30 pulsars could be radio- 
detected. A stud y of the dispersion measure of the observed pulsar pe- 
riod derivatives dFreire et alj|2001l) provided an estimation of the av- 
erage msPSR intrinsic period derivative with </'//'),„, ~ 10"'^ s"' and 
hence a surface dipole magnetic field B, » 2.6 x 10** G and a spin-down 
power of Ljrf a lO^'* er gs~'. 

At higher energies. lHeinke et alj d2005l) reported, from Chandra X- 
ray observatory data on 47 Tucanae, some 200 X-ray point sources, 
which belong to several object classes including cataclysmic variables, 
low-mass X-ray binaries (XRB), and the radio-detected msPSRs. They 
derive, from a tentative identification of the unknown sources they de- 
tected, an upper limit on the number of pulsars in the core of 47 Tucanae 
of about 60, assuming individual fluxes similar to the X-ray detected 
ones. Roughly two thirds of these msPSRs have a stellar companion 
(M < 0.2 Mq). The X-ray spectrum of a msPSR in a GC can be de- 
scribed by a thermal component plus single power law, with typical 
X-ray (0.5 - 6keV) fluxes around 10" ergs"' (Bogdanov et al. 2006). 
A few msPSRs exhibit X-ray pulsations, although with pulsed fractions 
below 50% for most of them (iCameron e t al. 2007). The presence in 
47 Tucanae of "hidden" msPSRs, detectable in hard X-rays but not 
in radio, has been excluded, within the uncertainty of the model by 
iTavanil ( |1991|) . by the high -energy X-ray (0.75 to 30MeV) upper limits 
reporte d bv COMPTEL ( fO'Flahertv et al.lll995h . From EGRET obser- 
vations, iMichelson et al.l ( Il994h produced a photon flux upper limit of 
5 X 10"^ cm"^ s"' above 100 MeV at 95% confidence level. In the same 
energy band, a detection of 47Tuc was just announced with the release 
of the FGST bright sourc e list^, see ^Abdo et al. (2009!), slightly below 
the EGRET upper limit bv lMichelson et al. ( 1994). These are discussed 
in Section[3] 

In the TeV range, previous observations of globular clusters re- 
sulted in upper limits. A limit on the steady photon flux fro m Ml 3 (5 
msPSRs, 7 kpc) was established by the Whipple Telescope dHall et al] 
120031) at 1.08 x lO^^cm-^s"' above 500 GeV. 47 Tucanae was ob- 
served by the Durham Mark III telescopes, with a resulting upper 
limit on the photon flux in pulsed emission from selec ted pulsars of 
4.4 X 10""cm"-s"', above a threshold of 450 GeV dBowden et alj 
Il991h . Periodic VHE emission from an XRB in 47 Tucanae above 5 TeV 
was reported once, by de Jager et al. (198^) during a remarkably high 
X-ray flux episode tAuriere et al, 1989) . Such event has not been re- 
ported since then in 47 Tucanae. 

The large number of identified msPSRs in 47 Tucanae and the com- 
pactness of the msPSR population at relatively close distance motivated 
H.E.S.S. observations of this GC, to investigate the predicted VHE 
emission from this class of objects. The results of these observations are 
presented in Section [2] Given the unknowns regarding VHE-emitting 
XRB, the interpretation given in Section[3]centers on a collective signa- 
ture from the msPSR population at TeV energies. 

2. Observations and Analysis 

H.E.S.S. is an array of four Imaging Atmospheric Cherenkov 
Telescopes, located in the Khomas Highland of Namibia. Stereoscopic 
analysis methods allow efficient background (cosmic ray) rejection and 

http://fermi.gsfc.nasa.gov/ssc/data/access/lat/bright_src_list/ 
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Fig. 1. Upper limit integral flux curve derived from the H.E.S.S. 
observations of 47 Tucanae (assuming a photon index of a = 
2), for "standard" cuts, at the 99% confidence level. Predicted 
fluxes for 100 msPSRs were added for com parison, rescaled for 
a dista nce of 4 kpc. (*) Curve adapted from lBednarek & Sitarekl 
(120071) . for ee = 0.01, = 100 GeV and a = 2, rescaled to 
Lsd - 10^"^ ergs"' (see Section|3]for details). 



accurate energy and arrival direction reconstruction for -y-rays in the 
range 100 GeV - 100 TeV. For point-like sources, the system has a de- 
tection sensitivity of 1% of the flux level of the Crab Nebula above 
1 TeV with a significance of 5 cr in 25 hours of observation. A thor- 
ough discussion of the H.E.S. S. standard analysis and performance of 
the instrument can be found in lAharonian et al] d2006ah . 

A total of 13 hours of 4-telescope data have been taken by 
H.E.S.S. between October and November 2005 on 47 Tucanae (ex- 
cluding data taken during bad weather or affected by hardware irreg- 
ularities). The target was observed with an average zenith angle of 
50 ° and mean target offset of 1 ° from the centre of the field of view. 
Ap plying the H.E.S.S. anal ysis "standard" cuts for point-like sources 
(see lAharonian et al.l2006al) . the energy threshold is about 800 GeV and 
the point-spread function above 1 TeV is 0.11°, too large to resolve the 
core of 47 Tucanae. Tighter cuts would slightly improve the sensitiv- 
ity and angular resolution but also increase the energy threshold, fur- 
ther reducing the chances of a detection according to the models (see 
Section[3]l. Several methods for y-ray reconstruction (the Hillas param- 
eters method and a semi-analytical approach described in lRolland et al.l 
l2004ft and background estimation (the "ring" and "reflected" algorithms 
discussed in Aharonian et al. 2006a) were used, with consistent results. 

We find no significant y-ray event excess over the estimated back- 
ground from the direction of 47 Tucanae. With standard cuts and using 
the "reflected" background estimation method, the significance of the 
excess in the 0.11° radius integration area is 0.7 cr. This allows us to 
set an upper limit on th e flux from the targe t regio n. We determined 
upper limits according to lFeldman & CousinsI d 1 9981) with a 99% confi- 
dence level, assuming a point-like source and a power law photon flux 
energy spectrum of index a = 2. The integral flux upper limit discussed 
here and shown in Fig. [T] was derived using the standard Hillas anal- 
ysis, consistent within 20% with cross-check analyses. Increasing the 
photon index to a = 3 does not modify the result by more than 20%. 
The upper limit on the integral photon flux in the H.E.S.S. energy range 
for this data set (800 GeV - 48.6 TeV, from the energy range of the col- 
lected events) is 6.7 x 10"" cm"^ s"' or ~ 2 % of the Crab flux. This 
translates into a limit on the energy flux in the same energy range of 
6.8 X 10^' erg s"' when placing 47 Tucanae at 4 kpc distance. We also 
investigated an extended region (0.2° radius), without finding a signif- 
icant excess. We do not discuss the extended case further due to the 
compact distribution of the msPSRs in 47 Tucanae and the generally 
weaker limits derived for extended regions. 
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3. Discussion 

The H.E.S.S. upper limit on the y-ray flux emitted by 47 Tucanae can 
be confronted with scenarios of VHE y-ray emission by msPSRs in- 
volving accelerated leptons in progressively larger regions: close to the 
pulsar, inside the pulsar wind nebula (PWN), at the boundary of the 
eventual PWN, or further away in the GC where pulsar winds may 
interact. The comparison to PWNs detected in the VHE range is also 
discussed. We only consider here average properties of the msPSRs in 
47 Tucanae, as summarized in Section [T] for populations of 23 (de- 
tected) or 100 sources. While observational results favor smaller num- 
bers, r esults from dynamical models of GCs, e.g. from llvanova et al] 
( |2008|) , suggest possibly larger populations. Unless stated otherwise, 
the following constraints scale linearly with the number of pulsars. 

The production of y-rays in the pulsar magnetosphere has been 
proposed in (at least) two different general scenarios, which consider 
different production sites: th e "outer gap" or the "po lar cap". In the 
"outer gap" model (see e.g. IChen & Rudermanlll993h , low values of 
the surface magnetic field (estimated from the spin-down rate) and pul- 
sar period, which define the conditions near the light cylinder, are be- 
lieved to generally prevent VHE emission fr om msPSRs. Althoug h the 
"polar cap" model (discussed for instance in lHarding et alj|2005l) does 
not have such restriction on the conditions for VHE emission, both 
classes of model predict the flux to drop o ff sharply between 1 and 
100 Ge V, as discus s ed for a single msPSR in iChiang & Romanil ( 1 19921) 
and in IWang et alj (12005 ) for a large population. The upper limit by 
EGRET jMichelson et alJI 19941) does constrain some of these models. 
Pulsed emission is also predicted, e.g. in I Venter & de Jagg 1 ll2008h for 
47 Tu canae, to drop before 100 GeV, below the limit bv lBowden et alj 
( Il99 ih . The Fermi detection will undoubtedly renew the discussion on 
these processes, but interpretation in the 20 MeV - 300 GeV band will 
be challenging, between potentially pulsed emission from one or more 
msPSR, confused or unresolved sources, and the overall steady emis- 
sion component, which could be tied to scenarios also valid at energies 
above our quoted threshold. 

The IC component, produced either in the magnetosphere or fur- 
ther away from the compact object, does extend to the e nergies con- 
sidere d here, but in most cases with only very low fluxes jBulik et aU 
20001). Still, whe n cons i dering populations of sources, as done in 
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Venter & de Jagej ( 120081) : I Venter etal.l | |2009|) , IC emission from the 
pulsar nebulae might reach observable levels. Their Monte Carlo simu- 
lations of msPSR populations, accounting for the observed range of pa- 
rameters {P, Pj,,,, viewing geometry), predict the cumulative flux from 
100 msPSRs, as illustrated in Fig.[T] The efficiency of the IC emission 
in the pulsar nebula increases with the strength of the nebular mag- 
netic field B (in relation with the increased confinement time) until 
losses by synchrotron radiation become dominant. For a given value 
of B, the H.E.S.S. upper limit can be normalized by the predicted flux 
per pulsar to obtain the maximum allowed number of msPSRs, as done 
in Fig. [2] Large msPSR populations are thus excluded, down to 80 ob- 
jects for S = 12 jiG. In this model, the prediction falls short of providing 
constraints for only 23 msPSRs. Using the limit on the magnetic field 
strength in the nebula - po st-shock - of the millisecond PSR J0437- 
4715 bv lZavlinetaP ( l2002h of S < ISfjC and possibly lower (see the 
discussion in that reference), and assuming similar properties for 100 
msPSRs, this would suggest B < 5 jiG in the average pulsar nebula. 

Another scenario for producing VHE y-ray emission relies on par- 
ticle acceleration at the shock discontinuity of a PWN. Thorou gh dis- 
cussions on pulsar winds can be found in iKaspi et al.l ( l2006h . In X- 
rays, pulsar wind emission from msP SRs has been observ ed, in the 
so-called "black widow" discussed in IStappers et aU j2003l) . with lu- 
minosities similar to those of canonical pulsars , but not in a GC. I n 
this object, as well as for the "Mouse" pulsar dGaensler et alj[2004h . 
there are indications of interaction with the interstellar medium, sug- 
gesting a bow shock geometry primarily driven by the pro per motion of 
the pu lsar rather than by its accelerated particles. However, I Cheng et alj 
( 120061) established that in a GC such bow shock emission would be 
hampered by the geometry and stellar density. VHE y-ray emission 
from several PWNs has already been detected and lCheng et alj ( |1986|) 
suggested that msPSRs host the same leptonic emission processes as 
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Fig. 2. Upper limit on the number of msPSRs in 47 Tucanae for 
a given a verage m agnetic field in the pulsar nebula, using the 
model by I Venter et al. (2009.) and the H.E.S.S. flux upper limit. 
The dashed lines indicate the number of observed msPSRs (23) 
and the 100 msPSRs hypothesis discussed here. 
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Table 1. Upper limits on conversion efficiencies from spin-down 
power. 



young pulsars like Vela X (290 pc, Lj ^ ^ 10^^ergs~'). W ithout assum- 
ing a particular emission process (see lHorns et alj|2006l for a hadronic 
VHE emission model for the Vela X PWN), we derive the flux expected 
if similar objects were located in 47 Tucanae. The VHE detection of 
the Vela X PWN (Aharonian et al. 2006b) gives an integral photon flux 



F{E > 800 GeV) ^ 1.5 x 10"" cm"- s"'. Scaling for the distance and 
spin-down power of the pulsar associated with the Vela X nebula to 
the pulsars in 47 Tucanae amounts to a factor 5.3 x 10"''. We cannot 
constrain this model, as 840 "Vela-like" msPSRs would be required to 
reach our flux upper limit. From pulsar properties and measured fluxes, 
it is usual to estimate the fraction o f the spin-down pow er converted to 
y-rays, 6,^, as compiled recently in lHessels et alj ( l2008h for the VHE- 
detected PWNs in the 1 - lOTeV energy band. For 47 Tucanae, this 
fraction is limited by N,, x e'j"' < 0.7. Any msPSR population with 
A',, > 23 gives e'^'" in the broad range of detected PWNs (8 x 10"'' to 
0.05). From this point of view, the msPSRs in 47 Tucanae cannot be 
distinguished from the much younger and more energetic pulsars de- 
tected through the VHE emission of their PWN. Detailed studies of the 

specificities of each PWN might clarify the picture. 

Nonetheless, a scenario bv iBednarek & Sitarekl ( l2007h proposes 
that the energy of primary particles for the IC process increases through 
the interaction of the leptonic pulsar winds inside the GC. No observa- 
tional evidence for such wind-wind interaction has yet been found. They 
predict appreciable VHE y-ray fluxes for a population of 100 msPSRs 
when the power emitted by each pulsar is fixed at 1 .2 x 10'^ erg s"' . The 
distribution in energy of the leptons produced by a pulsar is assumed to 
follow a power law of index a, above a minimum energy £,„,„ . In most 
cases, the predicted flux in the H.E.S.S. energy range for 47 Tucanae 
should be above the detection threshold of the instrument. According to 
this model, a non-detection translates in a limit on NpXe^, the number of 
pulsars times the conversion efficiency from the pulsar spin-down power 
into relativistic electron-positron pairs (and not e^j, from spin-down to 
photons). The available HESS d ata on 47 Tucanae do not all ow to reach 
the reference sensitivity used bv IBednarek & Sitarekl ( l2007h . estimated 
(for 50 h of observation at 20 ° zenith angle and 0.5 ° offset) as a photon 
flux of about 2.0 x lO"'"* cm"^ s"' above 800 GeV, a factor /,„,,, « 3.35 
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lower than the resu l t pres ented here. Besides, the values assumed in 
iBednarek & Sitare3 ( |2007|) for the distance to 47 Tucanae (4.5 kpc) and 
the individual spin-down power (1.2 x 10^^^ ergs"') may be too large. 

Overall, a factor /„,„., X (lBed„arem / ^d,u, ^ ^ ^jd„ty jBed,u,rek01 ^ 3J g 

must be applied when comparing their model predictions to the pre- 
sented H.E.S.S. upper limit. Since their original limit is on Np x 6^,, a 
linear rescaling can be applied when changing the number of pulsars. 
Rescaled conversion efficiencies, derived from the H.E.S.S. upper limit 
above 800 GeV assuming 100 msPSRs, are given in Tab. [T] for their 
model (noted :j:) and for the conversion from spin-down power to VHE 
emission (e'^'") discussed above. The comparison depends on the injec- 
tion spectrum of the leptons produced by the pulsars. All the proposed 
scenarios are constrained (e^ < 1) in the 100 msPSRs case, with most 
limits on the efficiency clearly below the estimated 6(, » 0.1 for the 
Crab nebula dBe dnarek & Sitare3l2007l) . even when assuming only 23 
msPSRs. The exception is the scenario where most of the leptons are 
produced with low energy (£,„;„ = 1 GeV and a = 3): the constraints 
weaken to 6^ < 0.6 for 100 msPSR and e,, > 1 (no constraint at all) for 
23 msPRSs. 



4. Conclusions 

The upper limit of the VHE y-ray photon flux obtained from H.E.S.S. 
observations of 47 Tucanae, F{E > 800 GeV) < 6.7 x 10"" cm^^ s"', 
is at present the second limit for a GC with a sizable population of 
msPSRs. Given the size of this population, it is the most constraining 
upper limit on the flux from an ensemble of msPSRs so far derived. 

Comparing this result to emission models, we considered msPSRs 
as the only potential y-ray sources in the GC. Owing to the high energy 
threshold of these observations, emission models for the pulsar polar 
region, generally predicting low fluxes at these energies, cannot be con- 
strained, except when assuming msPSR populations muc h larger than 
consid ered here (23 - 100). These numbers, according to I Venter et alj 
( I2009h . may however be sufficient for the total IC emission to reach flux 
levels where the number of pulsars can be limited, depending on the 
strength of the magnetic field in the pulsar nebula, down to 6 < 5 fiG in 
the average pulsar nebula for 100 msPSRs. The limit on the conversion 
efficiency from spin-down power to VHE flux (see Tab. [TJ is compat- 
ible with the results avail able for VHE-detected PW Ns. Collective IC 
emission as proposed by IBednarek & Sitare3 ( 120071 ) cannot be more 
efficient than in the Crab nebula for most of their sets of parameters. 
Complementary constraints at lower energy should follow the detection 
of 47 Tucanae by the Fermi Large Area Telescope, but given the possi- 
ble complexity of the emission in the GeV range, the connection to the 
VHE band cannot be assessed here. 
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